Ab initio calculations of the ^^S ^Pj and 
33g-^ 37,35qy 2po hyperfine structures. 

T Carette and M R Godefroid 

Chimie Quantique et Photophysique, CP160/09, Universite Libre de Bruxelles, 
Av. F.D. Roosevelt 50, B-1050 Brussels, Belgium 

E-mail: t caretteOulb . ac . be , mrgodef (§ulb . ac . be 

Abstract. We present highly correlated multi-configuration Hartree-Fock (MCHF) 
calculations of the hyperfine structure of the 3p^ ^Pj levels of ^^S~ and ^^'^'^Cl. We 
obtain a good agreement with observation. The hyperfine structure of the neutral 
sulfur "^^S 3p'^ ^Pj lowest multiplet that has never been measured to the knowledge of 
the authors, is also estimated theoretically. We discuss some interesting observations 
made on the description of the atomic core in MCHF theory. 

PACS numbers: 32.10.Hq,31.15.A-,32.10.Fn 

Submitted to: J. Phys. B: At. Mol. Opt. Phys. 

Keywords: negative ions, sulfur, chlorine, hyperfine structures 



1 February 2011 



Ab initio study of the ^^S ^Pj and ^'^S / ^^'^^C/ 3p^ ^Pj hyperfine structures. 2 

1. Introduction 

As the experimental techniques evolve, the window of isotope effects in negative ions is 
becoming wider and wider [Il|2l[3llll[5l|6|,[7]. It offers new possibilities for understanding 
negative ions, correlation effects and the effect of nuclear properties on the electronic 
structure. In order to keep to this agenda, we need accurate models. 

The hyperfine interaction, as arising mainly from the non-spherical nature of the 
electronic wave function near the nucleus, is an operator that is particularly sensitive 
to both core and valence correlation effects. It is also sensitive to contributions that do 
not affect the energy much. When trying to calculate it, we face the difficulty of getting 
simultaneously a good description of the valence correlation - crucial for negative ions -, 
core- valence and core correlation. It is therefore a useful tool for testing new approaches. 

Experimental studies of the hyperfine interaction in atomic negative ions are 
scarce. A pioneer study was made by Mader and Novick on the ^He~ Is2s2p ^P° 
metastable state |[8]. Much more recently, the hyperfine structure of the only known El 
transition between bound states of a negative ion, the Os~ , was studied [9] . In between, 
Trainham et al |T] performed a study of the ^^S~ ^-^3/2 hyperfine structure. The S~ 
ion offers a good possibility for testing the methods on a middleweight system. 

Recently, the isotope shift on the electron affinity of sulfur between the zero-spin 
isotopes 32 and 34 was measured and calculated 0. A similar study involving the isotope 
33 of sulfur, of spin / = 3/2, would require the knowledge of the hyperfine structure of 
the involved states [S]. However, if the ^^S~ ground state hyperfine structure is known, 
the one of the ^^S 3p^ ^Pj multiplet is not, to the knowledge of the authors. 

We present in Section |3] large scale MCHF calculations on the iso-electronic systems 
S~ and CI, testing the possibility of using the valence correlation orbitals [lOJ for 
describing core effects. This choice of orbital set was successfully adopted for including 
a core- valence contribution in the calculation of the isotope shift on the electron affinity 
of sulfur [6] and chlorine [7]. In Section |4| we calculate the hyperfine structure of 
the neutral sulfur lowest multiplet and point out a remarkable interplay of core and 
valence correlation in a MCHF description. We compare our theoretical results with 
experiment in Section [5] and conclude in Section [6] that some computational strategies 
provide a better definition for the core region. 

2. Theory 

2.1. The MCHF expansion 

The multiconfiguration Hartree-Fock (MCHF) variational approach is based on the 
following ansatz 



where the ^{'yiLSMiMs) are configuration state functions (CSF) [Tlj that are built 
on a single orthonormal basis set. In this theory, all Cj are optimized and some - 




(1) 
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if not all - orbital radial functions are varied. We use the ATSP2K implementation 
of this method [12]. In this work, we first perform MCHF calculations to determine 
the orbital basis set. We then do subsequent larger configuration interaction (CI) 
calculations - i.e. optimizing only the q. An orbital active set is defined as the set 
of all orbitals characterized by quantum numbers with n < Umax and / < Imax, and is 

denoted [''T-max^maxl • 

Our scheme is based on the concept of "multi-reference-interacting" (MR-I) set of 
configuration state functions detailed elsewhere [6]. It is defined as the set of all CSFs 
that interact to first order with at least one component of a given multi-reference. The 
LS angular momenta of the occupied sub-shells are coupled by decreasing n and /. 



2.2. Hyperfine interaction 

The hyperfine structure of a LSJ level is caused by the interaction of the electrons 
with the atomic nucleus of angular momentum quantum number /. The total atomic 
angular momentum quantum number is denoted F. To first order, the hyperfine energy 
correction has the form 

^) - ^'2 + 8/(2/-l)J(2J-l) ■ 

where Aj and Bj are respectively the magnetic dipole (Ml) and electric quadrupole (E2) 
hyperfine constants and C = F{F -|- 1) — J(J + 1) — /(/ + 1). The theory underlying 
the computation of hyperfine structures using MCHF wave functions can be found in 
references [131 Ull IE] ■ It is possible to express the non relativistic hyperfine interaction 
in terms of the J-independent orbital (a/), spin-dipole (a^), contact (oc) and electric 
quadrupole (6) electronic hyperfine parameters [I3] 

N 

ai 

i=l 
N 



(TLSMLMslJ^^o^i^y^'l^LSMLMs) , (3) 

i=l 
N 

(FL^Mz^M^I ^2Cj')(2)4^)(?)r-3|rL5MiM5) , (4) 

i=l 
N 

(FL^M^M^I ^24')(^)rr25(r,)|rL^MiM5) , (5) 

i=l 
N 

b = (rL5MLM5|^2Cj'^(z)r-3|rL5MiM5) , (6) 



ad 



an = 



i=l 

and calculated for the magnetic component Ml = L and Ms = S ^ 



2.3. The distinction between core and valence electrons 

When studying an atomic system, it is common to distinguish core and valence electrons 
based on the mono-configuration approximation. This separation can usually be made 
quite efficiently by choosing a rare-gas-like core. In our case, the core is neon-like 



Ab initio study of the ^^S ^Pj and ^'^S / ^^'^^C/ 3p^ ^Pj hyperfine structures. 4 

(ls^2s^2p^) and the valence is composed of the n = 3 electrons (3s^3p^ for S~ and CI, 
and 3s^3p^ for S). From that distinction, a lot of concepts based on a first-order picture 
of the correlation arise: core polarization, valence, core- valence and core-core correlation 
(V, CV, CC), etc. This terminology remains vague and the contributions those terms 
refer to often depend on the method pil| [TB| \T7\ ITS] . In MCHF theory it is in general 
impossible to have a clean partition of those different contributions due to systematic 
rearrangements of the {q} MCHF eigenvector and orbital shapes from one model to 
another. The variational procedure is not perturbative and the core-valence distinction 
is based on a perturbative picture of the correlation. Still, the core-valence separation 
has a physical background so that it should not be discarded. Hence, we need a reliable 
definition of the core in order to interpret the so-called "core effects" . 

In this work, the uncorrelated core is defined as being the A^c-electron core state 
built on the Hartree-Fock orbital set optimized for the total system. In case of multiple 
solutions in the HF model, we select the one satisfying the Koopmans requirement [19] . 

3. Hyperfine Structures of the "^Pj levels of S and CI 

3.1. Definition of the orbital set: MCHF calculations 

To explore the impact of the choice of the MCHF model defining the orbital set on the 
results of the open-core CI models, we use two references: the single main configuration, 
denoted SR, and the set of all valence single and double (SD)-excitations in the = 3 
shell, denoted MR, i.e. 

MR= {l,2}i°{3s,3p}^{3}l (7) 

We then perform HF frozen-core MCHF calculations on the SR-I and MR-I sets obtained 
by activating all the electrons in orbital basis sets ranging from [4/] to [13/i] and 
[lO/i] for SR and MR, respectively. The so-obtained one-electron radial functions are 
denoted SR-I-C and MR-I-C. In addition, we perform valence MCHF calculations on 
the MR-I("nZ], n = 4 — 13, / < 5 (/i orbitals). The resulting basis set is denoted MR-I-V. 
In all calculations, the [Ne] core is frozen to its HF shape. 

Tables [l] and [2] display, respectively for CI and S~, the weights of the first 
configurations in the various MCHF calculations with Umax^max = lO/i. The weight 
of a configuration is defined as 




where the sum runs over the CSFs belonging to the configuration. We observe that the 
order of the few most important configurations is very similar in all cases. Still, some 
differences appear mainly due to the fact that some radial functions Pni{f) become inner 
orbitals in the MR-I-C and SR-I-C models. This is illustrated in the first four columns 
of Table [3] in which we compare the mean radius of the spectroscopic and correlation 
orbitals of the MR-I-C model with the ones of the MR-I-V model for both CI and S~. 



Ab initio study of the ^^S Sp"^ ^Pj and ^^S / ^^'^^Cl 3p^ ^Pj hyperfine structures. 5 

Table 1. Sorted weights of the first configurations in the MR-I-V, MR-I-C and SR-I-C 

wave functions of CI. The active set is \10h] . The Is and 2s sub-shells are closed in 
all those configurations. 
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MR-I-V 




config. 
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config. 
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config. 
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0.1145 


2p63s23p34p2 


0.0649 


2/35^3^34^2 
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0.0527 
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^jJ kjO Ol* i / 


0538 


3s43T)43^i4fi 
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0.0691 


2p63si3/3rfi4pi 


0.0411 


2/339^3^2 


0.0416 


3si3p44si4pi 


0.0584 




0.0403 


2/3si3p'*3di4pi 


0.0407 


3s43p43(ii4pi 


0.0537 


2p^3p^3(f 


0.0346 


2p'^3s^3p^As^Ap^ 


0.0407 


3s23p44/i 


0.0475 


2p63s23p44ji 


0.0311 


2/3s23p44ji 


0.0323 


3p53ci2 


0.0417 


2p^3s^3p'^4p^5s^ 


0.0310 
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33233934/2 


0.0370 


2p^3s'^3p^3d^5g'^ 


0.0269 


2p63s23p33rfi5gi 


0.0273 


35233344391 


0.0322 


2p43s23p54^i5^i 


0.0268 


2^43523^54^15^1 


0.0267 


3s23p33d45c/i 


0.0306 


2^43523^55^2 


0.0257 


2^43523^55^2 


0.0253 


35233945^1 


0.0292 


2p63s23p34j2 


0.0252 


23363523^34^2 


0.0248 


3si3p44pi4di 


0.0271 


2p'^3s^3p^5p'^ 


0.0249 


2^43523^55^2 


0.0245 


3s23p34si4(ii 


0.0257 



Table 2. Sorted weights of the first configurations in the SR-I-C, MR-I-C and MR-I-V 

wave functions of S~. The active set is [lO/i]. The Is and 2s sub-shells are closed in 
all those configurations. 



SR-I-C 




MR-I-C 




MR-I-V 
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2^63523^34^2 


0.0848 


3523^34^2 
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0.0557 


2^63513^43^14^1 
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3s!3p43di4pi 


0.0755 
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0.0496 
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0.0500 
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0.0703 
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0.0441 
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0.0624 
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0.0424 
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0.0436 
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0.0606 


23563/3^2 


0.0342 
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0.0420 


3523p44/i 


0.0494 


235635I3354435I55I 


0.0328 
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0.0328 


3523p45pi 


0.0463 


23543523^55^2 


0.0309 
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0.0305 


3p53d2 


0.0418 
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0.0296 
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0.0295 


3523p34/2 


0.0355 
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0.0276 


2^63523^33(^15^1 


0.0281 


3s23p34si4di 


0.0343 


2^43523^55^2 


0.0263 


2^63523^45/1 


0.0272 


35i3p44pi4di 


0.0342 


2^63523^45^1 


0.0260 


2^43523^55^2 


0.0259 


35i3p43di5/i 


0.0335 



In open-core calculations, most of the n — h orbitals become inner orbitals and Pid-, Pif 
contract significantly. The SR-I-C orbitals do not differ strongly from the MR-I-C ones 
and are therefore not presented here. 



Ab initio study of the ^^S ^Pj and ^'^S / ^^'^^C/ 3p^ ^Pj hyperfine structures. 6 



Table 3. Mean radius (in oq) of the n = 1 to 5 orbitals of CI, and S in the 
MR-I-V[10/i] (valence) and MR-I-C[10/i] (open-core) models. Is, 2s, 2p are Hartree- 
Fock orbitals. 
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3p 
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1.82778 


2.28135 


2.29572 


2.02884 


2.03898 


3d 


1.76250 


1.75997 


2.15941 


2.16075 


1.95878 


1.96343 


4s 


2.03477 


1.92840 


2.44030 


2.39304 


2.40713 


0.71754 


4p 


2.41053 


2.13745 


3.34104 


2.80812 


2.81064 


0.70561 


M 


2.23120 


1.35932 


2.86818 


1.64193 


2.45573 


1.42118 


4/ 


1.82616 


1.30900 


3.03237 


1.17530 


2.98537 


1.20579 


5s 


2.14697 


0.69456 


2.71559 


0.74676 


2.28372 


2.13212 


5p 


2.16166 


0.68502 


2.88173 


0.75013 


2.42314 


2.38464 


5d 


2.08468 


0.54938 


2.04493 


0.76611 


1.70728 


0.60688 


5/ 


2.06913 


0.65051 


2.31986 


1.55795 


2.31174 


0.63876 


5g 


1.85751 


1.78797 


3.40818 


2.12699 


3.83176 


1.94307 



3.2. Fully correlated CI calculations 

The results of Tables [T] and |2] indicate that, even if the variational contents of open- 
and closed-core calculations are a priori very different [20], it is here possible to find 
extensions of the models developed above that remain comparable. Indeed, if we choose 
as multi-references the five or nine first configurations in the sorted lists 

MRs = {3s23p^ 3s23p^3rf^ 3s^3p^M\ 35^3^^^, ^s^^p^UHf] (9) 
MRg = MRs U {3s^3/4s^4p\ 3s^3p^3rf^4p\ 3/3c/^ ^s^^p'^Af} (10) 

for both systems, the selected MRs account for approximately the same correlation 
effects in all S~ and CI orbital sets. It allows a significant comparison between the 
CI calculations performed on the MR5- and MR9-I-C multi-reference-interacting sets 
performed with the orbital sets arising from either the SR-I-C, MR-I-C or MR-I-V 
MCHF calculations. 

Figure [1] and the two top plots of Figure [2] present the convergence of the a;, a^, 
Oc and h hyperfine parameters (|3]j6]) calculated with the MR5-I-C model in the three 
orbital basis sets (SR-I-C, MR-I-C and MR-I-V). All parameters are given in atomic 
units (cto^). In all cases the MR-I-V [13/;,] and SR-I-C [13/;,] results differ by about 
< 2 10^^ Oq^, the final Aj constants themselves differing by less than half a percent. 
The convergence of the calculations with the orbital sets obtained from the MR-I-V 
model is slow compared to the one of the calculations based on the SR- and MR-I- 
C orbital sets. It is remarkable to note that, to the contrary of the magnetic dipole 
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Table 4. Best estimates of the energy (in Eh) and hyperfine parameters (in Uq^) of CI 
(see text). The final values are the results of the MR9-I-C[12ft,] in which the impact 
of the 13^^ correlation shell on the MR5-I-C model is reported. 





MR5 




MRg 
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12h 
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460.12145 


-460.12315 




ai 


7.15768 
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7.14014 
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a-d 


-1.50471 


-1.50541 


-1.50369 


-1.50439 




0.49079 


0.48733 


0.52870 


0.52525 


b 


2.84768 


2.84738 


2.83912 


2.83882 



parameters (a), the electric quadrupole parameter (b) converges as rapidly in the SR-I- 
C and MR-I-V basis sets. 

Comparing the S~ and CI trends, it is the similarities that first strike. However, 
the j/-axis scales being the same between the left (S~) and corresponding right (CI) 
plots, we observe easily that (i) the S~ hyperfine parameters are shifted to smaller values 
compared to the CI parameters, and (ii) the S~ hyperfine structure calculations converge 
slightly faster that the CI ones. These two observations can be easily understood by 
the mere diffuseness of the negative ion electron charge distribution. Indeed, when 
an electron attaches to a neutral atom, the valence shells spread and the core-valence 
separation becomes larger. Then, even though the negative ion have more electrons, 
the core-valence overlap is only slightly larger in the negative ion than in the neutral, 
even if there are more valence electrons in the former than in the latter [71 [21]. Hence, 
the negative ion core is more spherical and the hyperfine constants tend to be smaller 
than in comparable systems. Consequently, even if negative ions are highly correlated 
systems, they are characterized by a comparatively smaller core-valence correlation. In 
particular, the hyperfine structure of negative ions is shghtly less sensitive to correlation 
effects than one would first expect. 

Even though the MR5 and MRg sets do not account for exactly the same correlation 
effects depending on the active set, the approach based on the orbitals obtained at the 
SR-I-C[13/i] level are superior. Tables [i] and [s] show the total energies and ai, a^, ttc and 
b parameters deduced from the MR5-, MR9-I-C[12/i] and MRs-I-CflS/i] models using 

th 

the SR-I-C orbitals. The final results are obtained by reporting the impact of the 13 
shell on the MR5-I-C model on the MRg based values. 

4. Hyperfine Structures of the 3p'^ ^Pj levels of S 
4.I. Definition of the orbital set: MCHF calculations 

For the 3p^ ^P state of neutral sulfur, we first select the references from which the 
MCHF expansions are generated. Similarly to CI and S~, we use the mono-reference 
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Figure 1. For both S~ (left plots) and CI (right plots), values of non-relativistic 
magnetic dipole hyperfine interaction parameters ai,ad and Uc obtained by CI 
calculations in the MR5-I-C model, in atomic imits (oq^); with the three explored 
orbital basis sets SR-I-C, MR-I-C and MR-I-V, as a function of the active space \nl], 
I < 5. 
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Figure 2. For both S~ (left plots) and CI (right plots), values of non-relativistic 

electric quadrupole hyperfine interaction parameter b (in a^^) and energy E (in E/j) 
obtained by CI calculations in the MR5-I-C model with the three explored orbital basis 
sets SR-I-C, MR-I-C and MR-I-V, as a function of the active space \nl], I < 5. 
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Table 5. Best estimates of the energy (in E^^) and hyperfine parameters (in ^) of S 
(see text). The final values are the results of the MRg-I-C[12ft,] in which the impact 
of the 13^^ correlation shell on the MR5-I-C model is reported. 
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Table 6. Sorted weights of the first configurations in the SR-I-C, MR-I-C and MR- 
I-V wave functions of S. The active set is [10ft]. The Is shell is closed in all those 
configurations. 
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0. 


.0339 


3523/4/1 


0, 


.0489 


2522/3323^24^15^1 


0, 


.0324 


2522/3513/3^15/1 


0. 


.0328 


3523/4^1 


0, 


.0486 


2522/3513/3^16/1 


0, 


.0312 


2522/3523/4/5^1 


0, 


.0324 


3513/3^15/1 


0, 


.0356 


2s22/3s23/4di5di 


0, 


.0289 


2522/3523/4^15^1 


0, 


.0285 


3523/5^1 


0, 


.0313 


2s22/3s23/4/i 


0, 


.0271 


2522^63523/4/1 


0, 


.0281 


3p54pi 


0, 


,0292 


2522/3523/5^2 


0, 


.0260 


2522/3523/5^2 


0, 


.0256 


3513^23^3 


0, 


,0282 


2522/3523p''4d2 


0, 


.0260 


2522/3523/4^2 


0, 


.0256 


3523/4/2 


0, 


,0252 


2512/3523/4514^1 


0, 


.0253 


2512/3523/4514^1 


0, 


.0250 


3si3/4]5i4di 


0, 


,0247 



(SR) and a multireference 

MR= {l,2}i°{3s,3p}^{3}l (11) 

The weights of the most important configurations in the SR-I-C, MR-I-C and MR-I-V 
calculations using the active set [lOh] are presented in Table [6] The corresponding 
mean radius of the orbitals are shown in the fifth and sixth columns of Table [H 
We observe that the eigenvector composition differences between the open-core and 
valence calculations are much more pronounced than in the S~ and CI cases. For 
instance, the ls^2s^2p^3s^3p"^4p^ configuration gains a significant weight in both open- 
core calculations, the Ap orbital being localized in the inner region of the atom. This 
important core effect is accompanied by a decrease of the closed-core configuration 
3s^3p^np^, n = 4 for MR-I-V and n = 5 for SR- and MR-I-C models. In the open-core 
models, the cumulative weight of these two configurations is 0.055, which is comparable 
to the weight of the 3s^3p^4p^ in the closed core model (= 0.068). It indicates that the 
correlation of the 3p^ and 2p^ electrons is similar so that, even if we fix the core to its 
HF shape, high-order correlation effects mix the core and valence electrons through the 
2p'^3p'^4p'^ configuration, {3p, 4p) being variational. Similar observations were made in 
the context of the calculation of the electron affinity of Chlorine |7|, in which case the 
negative ion neon-like core can strongly mix with the valence 3p® electrons compared to 
the neutral MCHF solution. 



Ab initio study of the ^^S ^Pj and ^'^S / ^^'^^C/ 3p^ ^Pj hyperfine structures. 11 

Table 7. Best estimates of the energy (in E/i) and hyperfine parameters (in a^'^) of 
the neutral sulfur state. The final results are obtained by reporting the impacts 
of the 13*"'^ correlation shell on the MR5-I-C model and the impacts of the extension 
of the MRg reference to the MRg reference using the active set \9h~\ , on the values of 
hyperfine parameters obtained in the MR6-I-C[12/i] calculation. 







MR5 




MRe 




MRg 


final 


9h 


12h 


13h 


9h 


12h 


9h 


E 


-398.07547 


-398.08254 


-398.08329 


-398.07655 - 


398.08363 


-398.07952 




ai 


5.13968 


5.12937 


5.12869 


5.14143 


5.13020 


5.13748 


5.12557 


ad 


1.08020 


1.07860 


1.07918 


1.08205 


1.08039 


1.08447 


1.08339 




0.74528 


0.85674 


0.85825 


0.80130 


0.91908 


0.82046 


0.93976 


b 


-2.05860 


-2.05035 


-2.04957 


-2.05942 


-2.05064 


-2.05571 


-2.04615 



4-2. Fully correlated CI calculations 

The large differences between the roles played by each correlation orbital in the open- 
core and valence calculations prevent us from comparing the two models as in the cases 
of S~ and CI. Although the results of the different open-core CI models are coherent, 
we limit our discussion to the SR- and MR-I-C results. 

Driven by the eigenvector composition analysis presented in Table |6} we choose the 
multi-reference 

MR5 = {3s23p^ 3s^3p^3d^ 3s^3p^3d\ 3sV^dHf\ 3p^3d^}. (12) 

Then, for evaluating the impact of the addition of a core-polarization configuration in 
the reference, we build the MR 

MRg = MR5 U {2s^2p^3s'^3pHs^4p^} , (13) 

and finally we define 

MRg = MRg U {2s^2/3s23p^V, 2s'^2p^3sV^d^^p^}- (14) 

We show in Table [t] the hyperfine parameters for the sulfur 3p^ ^P ground state using 
these MRp-I-C models (p = 5, 6, 8) for the largest possible active set and using the 
SR-I-C radial functions. The final results are obtained by reporting the impacts of 

th 

the 13 correlation shell on the MR5-I-C model and the impacts of the extension of 
the MRg reference to the MRg reference using the active set [9/i] , on the hyperfine 
parameter values obtained in the MR6-I-C[12/i] calculation. 

5. Comparison with experiment 

The three considered isotopes, ^^S, ^^Cl and ^^Cl, have a spin / = 3/2 and 
respectively a magnetic dipole moment of +0.6438212(14) fiB, +0.8218743(4) fj,B and 
+0.6841236(4) hb [22j. Their nuclear quadrupole moments are still best determined 
using theoretical inputs [23], as we will see below. 
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Table 8. Comparison of theoretical and experimental Aj hyperfine constants (in 
MHz) for the lowest multiplet of ^^S, ^^S~, ^^Cl and ^''Cl. Non-relativistic estimations 
(NR) computed from the final results of Tables [4| [S] and [T] Relativistic corrections 
estimated with a CI-RCI approach |2?1 |5S] . 

33s 33s- 35C1 37q 



M Ai A3/2 Ai/2 A3/2 Ai/2 A3/2 Ai/2 

Non-relativistic 133.58 0.33 92.48 483.30 168.52 846.09 202.45 1016.46 

+rel. corrections 130.09 -4.62 91.43 496.28 167.52 872.05 201.25 1113.22 

Experiment 91.49(9)" 170.69^ 205.05'' 

° Reference [1]. 
^ Reference l26l . 



Table 9. Nuclear electric quadrupole moments {Q, in barns) estimated with 
equation (IT6| for ^^S' , ^^C\ and ^^Cl. From those Q, we calculate the Bj constants 
of lower multiplet with our value of b{S ^P). Non-relativistic estimations (NR) 
computed from the final values of Tables |4j [5]and[7j Relativistic corrections estimated 
with a CI-RCI approach [21 [25]. 

33s 35C1 



Q B2(33S) Bi(33S) Q Q 



Non-relativistic 


-0.0655(6) 


-31.49 


15.74 


-0.082265 


-0.064833 


-f-rel. corrections 


-0.0657(6) 


-31.60 


15.80 


-0.081764 


-0.064438 


Sundholm and Olsen"^ 


-0.0678(13) 


-32.60 


16.30 


-0.08165(80) 


-0.06435(64) 



Reference [51 for ^^S and [17] for the chlorine isotopes 



The non-relativistic Aj constants computed using the final set of hyperfine 
parameters of Tables |4] , [5] and [7] are shown in Table [8j We estimate the relativistic 
corrections by running mono-reference non-relativistic and corresponding relativistic CI 
calculations using the Pauli approximation with the SR-I-C[9/i] orbital set [25]. 

The neutral sulfur Ai(^P) hyperfine constant is characterized by strong cancellation 
between the spin-dipole {Ad) and orbit (Ai) contributions. Indeed, in the MR6-I-C[12^] 
CI calculation performed with the SR-I-C orbitals, we find 

Ai = 105.04 MHz Ad = -110.74 MHz A^ = 6.28 MHz . (15) 

We realize from Table |4||5] and [T] that the contact term (Ac) is by far the less converged 
contribution, bringing the largest source of uncertainty (~ 1 MHz). It is unclear to 
which extent the ^^S~ theory-experiment excellent agreement is accidental. 

As far as the Bj constants are concerned, it is more relevant to tabulate the electric 
quadrupole moments (Q) obtained for S~ and CI from the formula [11] 
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where Gg = 234.96475 for obtaining Q in barns when b is in a^^ and Bj in MHz. The 
available experimental data are: 

3^/2^^^-) = 26.24(23) MHz [T] (17) 
^3/2(^^01) = 54.872 905(55) MHz j^S] (18) 
^3/2(^^01) = 43.245 245(55) MHz (19) 

We compare our Q values with previous works [2l[27j in Table[9j Since, to our knowledge, 
there is no measurement of the neutral sulfur hyperfine structure, we use our value for 
the b{S ^P) parameter for estimating the Bj constants of ^^S with each (5(^'^S) value. 

6. Conclusion 

We perform MCHF-CI and RCI calculations of the hyperfine constants of the 3p^ ^Pj 
multiplet of ^^S~ and ^s.stqj ^^^^ ^^le ^Pj multiplet of ^^S. We obtain good agreement 
with previous theoretical works [21 EZ] for the nuclear electric quadrupole moments of ^^S 
and '^'''^^Cl, and with the A{3p^ ^-^3/2) experimental values [HEH]. It appears that the 
contact contribution, the main source of uncertainty in our non-relativistic calculations, 
is ten times smaller in S~ than in CI. We interpret this as an effect of an increased 
separation of the core and valence regions in negative ions. 

We show that, for sufficiently large active sets, orbitals optimized in closed-core 
MCHF calculations reproduce the results of proper open-core MCHF calculations. This 
approach has a significant advantage: the core-valence distinction in frozen and closed- 
core MCHF calculations is much cleaner. It allows to minimize the high-order core and 
valence mixing and hence get a better comparison between calculations performed on 
different systems {e.g. S and S^). 
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